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The Sultana-Dyer solution of general relativity representing a black hole embedded in a special 
cosmological background is analysed. We find an expanding (weak) spacetime singularity instead 
of the reported conformal Killing horizon, which is covered by an expanding black hole apparent 
horizon (internal to a cosmological apparent horizon) for most of the history of the universe. This 
singularity was naked early on. The global structure of the solution is studied as well. 
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I. INTRODUCTION 

Two main areas of gravitational physics, cosmology 
and the study of black holes, come together in the at- 
tempt to find exact solutions of the Einstein equations de- 
scribing black holes which are asymptotically Friedmann- 
Lemaitre-Robertson- Walker (FLRW) instead of asymp- 
totically flat. With the exception of static solutions, such 
as the Schwarzschild-de Sitter black hole or its general- 
izations, these exact solutions describe dynamical black 
holes. What is meant exactly by "black hole" when the 
metric is non-stationary and the usual teleological event 
horizon is not present is a non-trivial question which has 
led to extensive research on dynamical horizons and their 
mechanics and thermodynamics (e.g., [3, H, H, 0] and 
references therein). The well-known black hole thermo- 
dynamics [5] hinges on Hawking's discovery of thermal 
radiation from black holes [6| , and this calculation relies 
on neglecting the backreaction of radiation on the space- 
time metric. The full treatment of backreaction and a full 
understanding of time-dependent dynamical horizons are 
still in the future. It interesting, therefore, to find and 
study exact solutions of the Einstein field equations de- 
scribing black holes, by which we mean central singulari- 
ties covered by an apparent horizon-these can be used as 
testbeds for various theoretical characterizations of sur- 
face gravity, energy fluxes, etc. Only a few such solutions 
are known. A complication arising when a black hole is 
somehow embedded in a cosmological background other 
than the de (anti-)Sitter one sourced by a cosmological 
constant, is that the cosmic fluid tends to accrete onto the 
central black hole. Forbidding this accretion flow results 
in a rather artificial condition and in the old McVittie 
solution Q. 

There are other motivations for pursuing exact cosmo- 
logical black hole spacetimes: one is the problem of the 
effect of the cosmological expansion on local dynamics 
(Q and references therein), which generated the McVit- 
tie solution Q ■ It has been realized @ that participation 
in the expansion of the cosmic substratum may be more 
difficult to achieve for strongly than for weakly bound lo- 



cal systems, hence one wants to look at the most strongly 
bound local system, the black hole. This approach has 
led to the solutions of Sultana and Dyer jlO( and those 
of [111 EH 03 ■ Recent interest in phantom dark energy 
and its properties (including thermodynamics) [3, Ha ] 
motivated the study of the backreaction due to the ac- 
cretion of cosmological phantom energy onto black holes 
and the possibility that phantom energy may violate Cos- 
mic Censorship [Ig|. Exact solutions describing black 
holes in a cosmological fluid may also lead to toy models 
for evaporating black holes. Finally, alternative gravita- 
tional theories such as f(R) gravity have received much 
attention recently as possible models of the present accel- 
eration of the universe without postulating exotic dark 
energy ([I3|, see [l8j] for a review and [l£| for shorter 
introductions to this subject). f(R) gravity and all the- 
ories of modified gravity introduced for this purpose pro- 
duce an effective time- varying cosmological constant and 
spherically symmetric solutions in these theories are not 
likely to be asymptotically flat but rather asymptotically 
FLRW. Black holes will not be Schwarzschild-like but 
dynamical: in fact, the Jebsen-Birkhoff theorem is not 
valid in these theories and spherically symmetric solu- 
tions will, in general, be non-stationary (see, for example, 
the dynamical solution proposed in [20| in the context of 
metric f(R) = R 1+s gravity). The interest in cosmologi- 
cal black holes is not confined to this class of alternative 
theories (which are, anyway, special cases of scalar-tensor 
gravity [3 Gil): interest has come from the possibility 
that inhomogeneities lead to local variations of the effec- 
tive gravitational constant in scalar-tensor cosmologies 
PH . Exact cosmological and time-dependent black holes 
are of interest also in higher-dimensional Gauss-Bonnet 
gravity 



gravity 



22 1 and arise from intersecting branes in super- 
23 . 
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In spite of all these motivations, only a handful of ex- 
act solutions describing dynamical black holes embedded 
in cosmological backgrounds are known, and even fewer 
are properly understood. Here we analyse the Sultana- 
Dyer solution [Io| which reserves a few suprises, and dis- 
cuss its singularities, apparent horizons, global structure, 
and physical interpretation. We adopt the notations of 
Ref. 
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II. THE SULTANA-DYER SOLUTION 

The Sultana and Dyer solution [10( is a metric of 
Petrov type D interpreted by its discoverers as a black 
hole embedded in a spatially flat FLRW universe with 
scale factor a(t) cx t 2 l z (where t is the comoving 
time). It is generated by conformally transforming the 



(S) 1 ! 

Schwarzschild metric g ab — > g a b = ^ 2 9 a V with the in- 
tention of of mapping the Schwarzschild timelike Killing 
field £ c into a conformal Killing field (for £ c V c f2 ^ 0), and 
the Schwarzschild event horizon into a conformal Killing 
horizon. The choice = a(t) = ao t 2 ^ (the scale fac- 
tor of a dust-filled k — FLRW universe) generates the 
metric in the form given by Sultana and Dyer (Toj 
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ds 2 



a 2 (77) 
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2m\ 2 Am 
— =- an H — — drjdr 
r r 



^ ] df 2 
r 



(2.1) 



where m is a constant, dSl 2 = d0 2 +sm 2 6d(p 2 is the metric 
on the unit 2-sphere, a(rf) = ao?? 2 = aot 2 ^ 3 . We use the 
metric signature opposite to that of [Toj and the symbol 77 
and t for conformal and comoving time, respectively, with 
dt = adrj (these notations are switched with respect to 
those of [l0| but follow standard practice in cosmology). 

The Sultana-Dyer metric (|2. 1|) is obtained by confor- 
mally transforming the Schwarzschild metric written in 
Painleve-Gullstrand coordinates. Of course, the line ele- 
ment (|2.ip can be cast in a form explicitly conformal to 
Schwarzschild written in the more familiar Schwarzschild 
coordinates. By introducing the new time coordinate i 
defined by 



dt = dt 



2ma df 



(2.2) 



the line element (12.11) transforms to 



ds 2 



2m 



dP 



a 2 df 2 



-11- 'f) ,i,r 



2m 1 

r 

df 2 



-y 2 m 



a 2 f 2 dil 2 



f 2 dn 2 



(2.3) 



which is manifestly conformal to Schwarzschild with con- 
formal factor a and dt = adfj. It is obvious that it reduces 
to a spatially flat FLRW metric as r — > +00. 
By using the isotropic radius r defined by 



(2.4) 



f = r ( 1 H 

V 2r 

the line clement (12.31) becomes 



ds 2 = 



1- 



(1 



dt 2 



2 dn 2 ) 

(2.5) 



in isotropic coordinates (i,r,6,<p). The metric (|2 .5[) is 
recognized to be formally the same as the McVittie solu- 
tion [7[ but with the important difference that the metric 
coefficient m is now a constant, contrary to the McVittie 
case in which dm/mdi — — da/adi. We will refer to this 
equation as the "McVittie no accretion condition" . The 
latter stems from the explicit requirement that there is no 
accretion of cosmic fluid onto the central object. In mod- 
ern language, this has been recognized as the constancy 
of the Hawking-Hayward quasi-local mass M = m(t)a(t) 
(24| . Since the McVittie no accretion condition is vio- 
lated by the Sultana-Dyer solution, the latter describes 
an accreting object, as Sultana and Dyer recognize from 
inspection of the field equations and the fact that Tq 7^ 0. 

A second important difference between the Sultana- 
Dyer and the McVittie metric is that the latter is sourced 
by a single perfect fluid while the material source of the 
former is a mixture of two non-interacting perfect fluids, 
a null dust and a massive dust [10(. The stress-energy 
tensor is 



-Lab — -L ab 



(II) 



ab 



(2.6) 



where — pu a Ub describes an ordinary massive dust 



and T, 



(ii) 



Pn k a kb describes a null dust with density 



p n and k c k c — [lOj. The use of two fluids instead of one 
follows from the fact that, once accretion is allowed and 
the energy density and pressure depend on the radial, as 
well as the time, coordinate, solutions sourced by a single 
perfect fluid do not exist, except for the Schwarzschild- 
de Sitter metric (which is a special case of the McVittie 
solution) [Tim. 

We will make use of the quantity 



M(t) = ma(t) 



(2.7) 



which is not constant in the Sultana-Dyer solution. These 
authors identify a conformal Killing horizon, the locus 
f = 2m (or r = m/2 in terms of the isotropic radius 
r) [Tfj| . This is obtained by mapping the Schwarzschild 
event horizon of the seed metric (a Killing horizon) into 
the 2-surface f = 2rn. Therefore, Sultana and Dyer in- 
terpret their metric as describing a black hole embedded 
in a special (spatially flat) FLRW background universe. 
The authors identify certain problems with this solution: 
the accretion flow onto the black hole becomes superlu- 
minal, and the energy density of the cosmic fluid becomes 
negative after a certain time near the conformal Killing 
horizon at f = 2m. They proceed to study radial null 
geodesies and surface gravity on this surface, and the be- 
haviour of timelike geodesies representing, e.g., planetary 
orbits, in this metric. This is inspired by the old prob- 
lem of general relativity of whether the cosmic expansion 
affects local systems (see [8( for a recent review). 

Here we take the study of the Sultana-Dyer solution a 
step further. We find that it is quite different from the 
picture provided by these authors: in particular, we show 
that the locus f = 2m is not a conformal Killing horizon 
but a spacetime singularity. It is a weak singularity in 
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the sense of Tipler's classification [25f. If the current 
interpretation of the McVittie metric [24} applies to the 
Sultana-Dyer solution, the two disconnected regions f > 
2m and < f < 2m correspond to a black hole and white 
hole region, respectively. The spacetime singularity f = 
2m is covered by an apparent horizon, which we locate 
together with a cosmological horizon. This alleviates the 
problems of negative energy density and superluminal 
flux near f = 2m. The singularity expands with time, 
comoving with the cosmic substratum but the horizon 
expands at a slightly smaller rate and, in the infinite 
future, it approaches the singularity asymptotically. The 
f > 2m Sultana-Dyer spacetime has also a Big Bang 
spacelike singularity. In the early universe, near the Big 
Bang, the f — 2m singularity was naked, and got covered 
by an apparent horizon only later. All these features 
resemble certain characteristics of the Fonarev solution 
[H or of generalized McVittie solutions [U, EI El • 



III. SINGULARITIES AND APPARENT 
HORIZONS OF THE SULTANA-DYER 
SOLUTION 

Let us begin by examining the locus f = 2m: this is 
not a conformal Killing horizon as stated in pjj |. but a 
spacetime singularity. In fact, the Ricci curvature is 



pa 



12 



6 (dfH + 2H 2 



a 2 (fj) 1 



2m \ 



1 



2m 



(3.1) 



(where H = d\na/dt) and it diverges as r — > m/2 (the 
discrepancy with the Sultana-Dyer paper is discussed in 
Appendix A). This is a covariant statement and not an 
artifact of the coordinate system adopted because the 
Ricci scalar is an invariant of the Riemann tensor. R a a 
also diverges as a — > 0, the usual Big Bang singularity of 
the background FLRW universe. The metric determinant 



-a (if) r 



- 4 sin 2 i 



(3.2) 



and does not tend to zero as f — > 2m, hence an object 
does not get crushed to zero volume as it approaches 
this singularity, which is a weak one in Tipler's classifi- 
cation 25]. It will be shown below that this singularity 
is covered by an apparent horizon and therefore, as in 
the McVittie metric, the regions f > 2m and < f < 2m 
may be interpreted as describing a black hole and a white 
hole region J24[. Although this interpretation does not 
seem absolutely compelling to us in view of the fact that a 
physical object or particle could potentially cross a weak 
singularity, we do not address this issue here and refer 
the reader to the comprehensive discussion of Nolan [24T ] 
and to the references therein. 

The issue is now to determine whether the r = 2m sin- 
gularity is naked or covered by an apparent horizon and, 
therefore, whether it really describes a black hole em- 
bedded in a cosmological background. Since the metric 



is non-stationary, there is no event horizon and the ap- 
propriate horizon notion is that of apparent horizon [lj. 
We proceed to rewrite the Sultana-Dyer solution in the 
Nolan gauge, in which it will be straightforward to de- 
cide whether apparent horizons exist and, if so, to locate 
them. 

Using the area radius R = af, which is a geometric 
quantity, eq. (|2.7p . and the fact that 



dR 

df = HRdr) 

a 



(3.3) 



where H = a/a and an overdot denotes differentiation 
with respect to t, the line element (|2.ip becomes 



ds 2 



2M H 2 R 2 \ , , 

1 - -n sir \dt 2 + - 

R 1-f 1 



dR 2 



2M 
R 



2HRdtdR 



1 



2M 



R 2 dn 2 



(3.4) 



in coordinates (t, i?, 9, ip), with the singularity located at 
R = 2M. Let us use A = 1 — 2M/R and the new time 
coordinate T defined by 



dT = — \ dt H — — ^-tt 

F V A 2 - H 2 R 2 



dR 



(3.5) 



where F(T(t, R), R) is an integrating factor satisfying 



d 



d 



dR \F J dt \F (A 2 — H 2 R 2 



HR 



(3.6) 



to guarantee that dT is an exact differential. After 
straightforward manipulations the line element (|3.4p is 
recast in the Nolan gauge as 



ds 2 



2M H 2 R 2 \ , , , 

1 TTTT F dT 



+ 



R 
dR 2 



1 - 



2M 



1 2M H 2 R 2 

1 - — - T^W 



R 2 dn 2 



(3.7) 



The apparent horizons, if they exist, are the locus g 
0, or 



(3.8) 



We discard the lower sign in eq. (|3.8p which corresponds 
to a contracting universe, and the apparent black hole 
and cosmological horizons are given by 



Rbh(t) = 



1 - VI - 8MH 
2H 



R c (t) 



1 + VI - SMH 
2H 



(3.9) 

respectively. It must be M < ^-g—, or d < for these 



apparent horizons to exist. In the limit m — > the black 



4 



hole and its apparent horizon disappear and the cosmo- 
logical apparent horizon reduces to the familiar surface 
R = -h. The fact that the proper radius of the latter 
is larger than R c can be interpreted as the pull of the 
central black hole on the cosmic fluid. 

Since R(l — HR) = 2M > at the apparent horizons 
and H > in an expanding universe, it is 



Rbh = 



2M 



1 — HRbh 



> 2M 



(3.10) 



and the singularity R = 2M(t) = 2ma(t), at which 
the flow is superluminal and the energy density becomes 
negative-definite, is hidden by the apparent horizon. This 
alleviates somehow the problems of this solution reported 
in However, since a(t) = dot 2 / 3 in comoving time, 
as t -> +oo, MH = ma = jffi -> and R bh -> 2M 
approaching the singularity asymptotically. 

By using the scale factor a(t) = aot 2 / 3 of the Sultana- 
Dyer metric it follows that the constant m must satisfy 
for the apparent horizons to exist. This condi- 



m < 



31 1 



16ao 



tion is violated at early times, implying that the R = 2M 
singularity was naked at early times and later on (at 

t = ("^f 1 ) 3 ) an apparent horizon appeared that immedi- 
ately bifurcated into a cosmological horizon and a black 
hole apparent horizon covering the singularity. 

Tha gobal structure of the Sultana-Dyer solution is 
analyzed in the next section. To conclude this section, 
we discuss the Misner-Sharp [28[ and Hawking-Hayward 
[29l |3t3 ] quasi-local energies. 

The Misner-Sharp mass Mms is defined in terms of 
the area radius R by 1 - 2Mp = ~V C RV C R [H which, 
at the black hole apparent horizon, yields 



Rbh 1 - VI - 8MH 
M MS = — = — 



(3.11) 



In order to compute the Hawking-Hayward quasi-local 
mass we introduce the null coordinates (u, v) defined by 



du = —= 
V2 




H 2 R 2 

1 2M 
1 R 



FdT 



dR 



1 



2M 
R 



H 2 R 2 



(3.12) 



dv = — = 



1 

V2 



2M H 2 R 2 



•l - 



dR 



R 



1 _ 2M 
1 R 



FdT 



1 _ Ml - 
1 R 1 



H 2 R 2 



(3.13) 



in terms of which the metric is reduced to Hayward's 
standard form [30( 



ds 2 



-2dudv + R 2 dn 2 . 



(3.14) 



Since dR 



1 



the Hawking- 



2M _ H 2 R 2 (dv-du) 

r t^m: V2 ' 
Hayward quasi-local mass is easily computed using the 
prescription for spherical symmetry [30l | 

M HH = R (RuRv + (3.15) 
which yields, at the black hole apparent horizon, 

M HH - Mm s - ^ (3.16) 



and coincides with the Misner-Sharp mass. Both masses 
diverge in the limit R — > 2M. 



IV. GLOBAL STRUCTURE 

The causal nature of the singularity and the black hole 
apparent horizon are determined as follows. The sin- 
gularity is characterized by the equation R — 2M = 
and the normal is obtained by taking the gradient of this 
equation (the limit R — > 2M is implicit in the following) . 
Unfortunately, the integrating factor F appearing in the 
line element (13.71) is not determined explicitly, hence it 
is more convenienent to use the coordinates (t, R, 9, ip) in 
which the metric is given by eq. (|3.4[) and its inverse by 



(9 ab ) 



-1 



-HR (-, 

1 R \ 



V 



-HR. 
2M 



2M 

R. 



H R 



1 

TP 







R 2 E 



-J 

(4.1) 

One obtains n a = V a i? — 2mdS a Q = (-2MH, 1, 0, 0) and 



g n b = 



-HR, -H 2 R 2 



2M 
1-jpO.O 



(4.2) 



and, in the limit R - 2M -> n a = (-2MH, 1, 0, 0), 
n a = (-2MiJ,-4M 2 i/ 2 ,0,0), while n c n c = 0. The 
R = 2M singularity is null, which is not surprising if 
one remembers that is is obtained by conformal mapping 
of the null event horizon of the Schwarzschild black hole 
used as the seed metric in this solution-generating tech- 
nique [Tol ] . 

The black hole apparent horizon has equation 
f(t, R) = 2RH - 1 + y/l - 8MH = 0. The normal has 
the direction of the gradient V a /, or 



(B(t,R),2H,0,0) 



where 



B(t, R) = 2RH 



4M [H 2 + H 



y/l - 8MH 



(4.3) 



(4.4) 
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while 



n a = g ab n b = - 



(B + 2H 2 R) 



2M 
R 



HRB ( 2M H 2 R 2 \ \ , % 

__ + 2 tf(l-_-___) !0 ,0) (4.5) 



1 



Therefore, 



2M 

Rbh 



2HRbh 



AM [H 2 +H 



VI - 8MH 



-AH 2 R 



2HRbh — 



AM if 2 +H 



VI - 8MH 



1 - 



1M 



(C 2 + AH 2 R) 



1 - 

\C\ 
HRm 



1 M 



(C + 2H 2 R,y -AH 4 R 2 



(C + AH 2 R) , 



(4.6) 



where 



R=2M 




Big Bang 



FIG. 1: The conformal diagram of the Sultana-Dyer space- 
time. The horizontal wiggly line at the bottom describes the 
Big Bang singularity, the wiggly line at 45 degrees denotes 
the R = 2M null singularity, and the solid line at 45 degrees 
describes future null infinity. Null geodesies end at future null 
infinity or at the black singularity (either when it is naked if 
started early on, or crossing the timelike black hole apparent 
horizon labelled AH). 



null singularity R — 2M (the conformal mapping of a 
Schwarzschild event horizon) which was naked at early 
times, until it got covered by a timelike apparent hori- 
zon appearing at a finite time. The conformal diagram 
is sketched in Fig. 1. 



C 



=£{l- VI - 8MH) - AMH 3 



if VI - 8MH 



(4.7) 



Now, 

C + AH 2 R = 



1 



if VI - 8MH 
+12MH 3 - 2H 2 ( 1 - Vl - 8MH) 



— ( I - \ I -Mff/ 



(4.8) 



and at late times if = J| and H = ^| , hence the numer- 
ator of eq. (|4.8[) is dominated by the term 12MH 3 and 
n a n a > 0: the black hole apparent horizon is timclikc at 
late times. At early times, but after the horizons have 
appeared, it is 



1 



if VI - 8MH 

-2H 2 (l - Vl - 8Aiff) 

12mcto 1 
~ < 7 / 3 if Vl - 8Ai# 



TT 

— ( 1 - Vl -8MH) " + 12Mif 3 



> 



(4.9) 



and the horizon is timelike early on. 

The global picture of the Sultana-Dyer spacetime that 
emerges from these considerations is that of a universe 
with a spacelike Big Bang singularity at t = 0, the 



V. DISCUSSION AND CONCLUSIONS 

There are many motivations, discussed in the Introduc- 
tion, to find and study exact solutions of the field equa- 
tions of general relativity or alternative gravity theories 
describing dynamical black holes embedded in a cosmo- 
logical background. Among these the Sultana-Dyer met- 
ric deserves some attention in view of its relative sim- 
plicity and of the technique used to generate it, which 
may lead to a wider class of solutions. There is evidence 
that among the class of generalized McVittie solutions 
[I2I H(| those conformal to the Schwarzschild spacetime 
("comoving solutions") are generic [HI], and this consti- 
tutes further motivation to understand solutions seeded 
by the Schwarzschild metric, such as the Sultana-Dyer 
one. 

To summarize and discuss the salient features of the 
Sultana-Dyer spacetime, we found that the conformal 
image of the Schwarzschild event horizon is not a con- 
formal Killing horizon but rather a spacetime singularity 
at R = 2m (where R is the area radius) . This singularity 
is weak in Tipler's sense and is null, as should be ex- 
pected because the seed Schwarzschild event horizon is a 
null surface. This expanding singularity is interpreted as 
the effect of the pull of the cosmological matter stretch- 
ing the r = singularity of the Schwarzschild spacetime 
into one of finit radial extent. Early on, this singularity is 
naked and only later a timelike apparent horizon develops 
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which immediately bifurcates into a black hole apparent 
horizon covering the singularity and an apparent cosmo- 
logical horizon. The radius of the cosmological horizon 
is smaller than the value if -1 of the Hubble radius in a 
FLRW spacetime, which may be interpreted as an effect 
due to the gravitational pull of the central object on the 
cosmic fluid. 

The singularity expands comoving with the cosmic 
substratum, while the apparent horizon expands at a 
slightly smaller rate and eventually comes to coincide 
with the singularity in the infinite future. Sultana and 
Dyer report superluminal flow near the singularity, and 
one may question whether the notion of black hole makes 
sense at all in the presence of superluminal flows which, 
in principle, allow particles to escape from the apparent 
horizon. However, it should be noted that this radial flow 
is always directed inward and nothing actually escapes 
from the black hole apparent horizon. In our opinion, 
this unpleasant superluminal feature is due to the sim- 
plicity of the model under study and will not be present 
in more realistic and sophisticated models yet to come. 
For the moment we content ourselves with understand- 
ing the simple model (|2.1|) . The Sultana-Dyer spacetime 
exhibits the spacelike Big Bang cosmological singularity, 
the null black hole singularity and, later on, the two ap- 
parent horizons. As is clear from the conformal diagram 
of Fig. 1, the timelike black hole apparent horizon meets 
the null singularity in the infinite future, with the dust- 
dominated universe expanding and diluting forever. 

In addition to the superluminal flow, which we do not 
regard as a serious flaw, the Sultana-Dyer solution ex- 
hibits other less desirable features: the cosmological fluid 
has negative energy density at late times near R = 2M 
(To| ; it would be desirable if the cosmological matter were 
described by a single fluid composed of particles following 
timclike geodesies instead of an odd mixture of a null and 
a massive dust; and limiting the scale factor of the uni- 
verse to the special choice a oc t 2 l z seems too restrictive. 
Work is in progress to find new exact solutions, both in 
general relativity and in alternative theories, which do 
not share these problems. 
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Appendix A 

Eq. (13. ip for the Ricci curvature shows a divergence 
at r = 2m, which is not noticed in the Sultana-Dyer pa- 
per. Here we discuss the likely cause of this fact. The 
Ricci scalar coincides (up to the constant K — 8ttG) 
with the negative trace of the energy-momentum ten- 



sor T ab and the matter energy density and (zero) pres- 
sure in the Sultana-Dyer paper are identified using their 
Einstein equations (13) with the stress-energy tensor 
given by their eq. (4). The latter is meant to pro- 
vide the stress-energy tensor T ab obtained by conformally 
transforming the vanishing stress-energy tensor T ab of 
the Schwarzschild solution, which we want to reconsider 
here. Under the conformal transformation of the met- 
ric g ab —> gab — ^ 2 9ab the Ricci tensor and Ricci scalar 
transform according to 

Rab = Rab - 2V„V 6 In n - g ab g ef V e V f In 

+2 (V a In SI) (V fc In SI) - 2g ab g ef (V e In SI) (V f In SI) , 

(A.l) 

K = [R\ - SE*1) , (A.2) 

respectively Q. Using the Einstein equations R ab — 
\g a bR = KT ab and R ab - \g a bR = nT ab , one obtains 

- 2V a V 6 ^ 4V a SlV b Sl 

g cd W c SlW d Sl n g cd V c V d Sl 
-g ab + 2g ab y - . (A.3) 

Sultana and Dyer instead have, for vanishing T ab (eq. (4) 
of0), 

W 2 Sl 2V„V h fJ 3 - 
KT ab = 2~g ab — a —^- - ~ 9ab ~g™V m SlV n Sl . 

The discrepancy between these formulas, plus the fact 
that the tilded operator V is not the correct one to be 
used in the expression of T ab , is likely to be the reason 
why the singularity at f = 2m is missed in [10]. The 
correct Ricci scalar can be calculated using eq. (|A.2[) with 
R a a = and a = i 2 / 3 = fj 2 yielding 

~ a msi en (fj 2 ) 

R a — fiT = ^g— = 2 9 6 0a 

12 - 12 (A5) 

ff(l~^) a 3 (l~^) • 1 ^ 

The singularity at f — 2m was also noted recently in [3ll ] . 

The Sultana-Dyer solution can be reobtained using as 
material source a single imperfect fluid with a radial 
(spacelike) energy current instead of a two-fluid mixture, 
as shown in the first of Refs. [ll|]. Eq. (90) of this work 
for G a a = —R a a is clearly singular at r = 2m. 
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